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Abstract 

The subject of low energy photon-photon scattering is considered in arbi- 
trary dimensional space-time and the interaction is widened to include scatter- 
ing events involving an arbitrary number of photons. The effective interaction 
Lagrangian for these processes in QED has been determined in a manifestly 
invariant form. This generalisation resolves the structure of the weak-field 
Euler-Heisenberg Lagrangian and indicates that the component invariant func- 
tions have coefficients related, not only to the space-time dimension, but also 
to the coefficients of the Bernoulli polynomial. 

1 Introduction 

That an interaction between electromagnetic fields and the vacuum of the Dirac field 
would allow such nonlinear processes as the scattering of photons by photons has been 
known for some time. Even before the emergence of the full scale theory of quantum 
electrodynamics (QED) Halpern, Euler, and Heisenberg, among others, discussed the 
phenomenon in the mid 1930s [1], [| 0, |J. However, it was in the early nineteen fifties 
that the subject was considered in detail using QED, and Karplus & Neuman [||. [| 
determined the full photon-photon scattering cross-section. 

The work presented in this paper is related to a desire to generalise the discussion 
of photon interactions to arbitrary dimensional space-time and to consider how the 
parameters of the interaction change in this case, and also with the number of scat- 
tering photons. This will be done by investigating the underlying effective Lagrangian 
for the interaction, and little consideration will be given to measureable cross-sections. 
The reason for this being that the cross-section determined by Karplus & Neuman, 
and others subsequently [0], for four photon scattering in four dimensional space-time 
is extremely small (i.e. a ~ 10 _30 cm 2 ). This is perhaps not unexpected for such an 
interaction, but it does mean that obtaining experimentally testable values is of lim- 
ited importance at this stage. The generalisation to arbitrary dimensional space-time 
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serves a twofold purpose. Firstly, it is useful to know the form of the result in the con- 
text of dimensional regularisation. More importantly perhaps it will provide greater 
insight into the structure of the Lagrangian in the dimensions of practical interest, 
i.e. D = 2, 3 and 4. As a scattering process involving the fewest number of photons 
(i.e. four) will have the largest cross-section, generalisation to interactions involving 
an arbitrary number of photons will also primarily be of interest in further resolving 
the structure in the Lagrangian. 

A scattering process involving N photons can be represented by a sum over permu- 
tations of an iV-photon electron loop. i.e. a coupling of N photons via an intermediate 
loop of electron and positron propagators. Direct evaluation of the Feynman ampli- 
tude in this case would be quite involved and this problem is most easily tackled 
by taking, as a starting point, the full 1-loop Lagrangian for the field involving a 
summation of all possible couplings between the electromagnetic field and a virtual 
matter field. For reasons of clarity it is useful to split the discussion of iV-photons 
interactions into two cases: those involving an even number of photons; and those 
involving an odd number. We shall proceed now to consider the first of these. 



2 2N-Photon Case 



Interactions involving an even number of photons particularly lend themselves to the 
term photon-photon scattering, as a situation involving the same number of initial 
and final photons can easily be envisaged. In the discussion which follows we shall 
implicitly assume that the interactions involve 2N photons, where N is integral, and 
thus the number of interaction photons involved is always even. In order to discuss 
these cases we can recall the 1-loop electromagnetic field Lagrangians obtained by 
Brown & Duff |§ for the cases of bosonic and fermionic virtual particle loops: 
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These results being generalisations of the four-dimensional Lagrangians obtained by 
Schwinger || 10]. Brown & Duff used momentum space functional methods to con- 
sider the Lagrangian in any field theory due to 1-loop quantum vacuum effects. The 
reader is directed to their paper for details of the derivation, while we shall simply 
use their results to consider the case of 2iV-photon scattering. It should be noted that 
these expressions are valid in strong fields as they contain all orders of the coupling 
to the electromagnetic field. Taking only the lowest order terms will correspond to 
the situation of weak fields. 

As it is a simple starting point we shall first consider the situation of photon 
scattering via a virtual meson loop (i.e. spin zero propagator particles). Techniques 
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used in this case will then be applied to the slightly more complex process in which 
photon scattering is mediated by an electron/positron loop. 



2.1 Scalar Case 

Our starting point here will be the 1-loop Lagrangian of Brown & Duff (Eq. |l]) which 
can be written in the following form 
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where 
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This Lagrangian, although it must be Lorentz covariant, is not explicitly so. Our task 
will be to extract a given order of the interaction (i.e. 2iV-photon scattering) and 
express this part of the Lagrangian in a manifestly covariant form. To this end we 
can follow the procedure described by Delbourgo & Matsuki [[II] and suppose that 
the function f(z) possesses a series expansion of the form 
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which is clearly true in the case above. Then writing the exponential and logarithmic 
functions as power series we can obtain 
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where new constants a m have been introduced to simplify the expression. As this 
relation is now in the form of a single power series expansion we can pick out any 
given term, and in particular the general term, which corresponds to a 2iV-photon 
scattering interaction, will be given by 
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In this expression the summation over rti runs from zero to a certain integer such that 
the partitioning constraint, 



rii + 2n2 + 3n3 + • • • + Nn 
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is satisfied. We should note here that implicit in this selection of a particular or- 
der of the expansion as a valid Lagrangian is any renormalisation of charge, and 
corresponding scale changes for the fields, which may be needed. 
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This Lagrangian is now in a manifestly invariant form and the problem has reduced 
to determining the coefficients a m . Fortunately a method for determining these factors 



has been developed in a rather different area of field theory. Delbourgo & Matsuki [1 I 
when considering gravitational and chiral anomalies [12|, |13[ obtained the following 
expression in terms of a two-form matrix version of the Riemann tensor R 
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As we see this is identical to the second exponential term in the 1-loop Lagrangian 
(Eq. ^) apart from the use of R in place of F. The fact that the same mathematics 
has turned up in apparently unrelated areas of field theory is certainly of note, but 
leaving this interesting point aside we can use the following expression, obtained by 
Delbourgo & Matsuki [lip, to determine the coefficients a m . 
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This general result can be applied to the particular case of the spin-0 Lagrangian 
in which we have 
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Determination of the coefficients can then be achieved by performing the integration 
and noting that 
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where B n are the Bernoulli numbers. We obtain 
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Thus, for a 2iV-photon process mediated by a virtual meson loop, the interaction 
Lagrangian is given by 
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We can now proceed to perform the integration, obtaining 
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The integral in this expression is simply the defining relation for the T-function, i.e. 
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We will find that it is particularly illuminating in the spinor case to use the the 
coefficients of the Bernoulli polynomial B n (x) [l 
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in place of the Bernoulli numbers and for later comparison we shall introduce them 
here. With the substitution B n (0) = B n the invariant coefficients become 

^ J^m+l22m 

aL = Am(2m)\ B2m{ ° 1 
and we can write the effective interaction Lagrangian in the following compact form 
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This is our final expression for the Lorentz invariant low energy effective inter- 
action Lagrangian for the electromagnetic field due to spin-0 quantum loop effects. 
Although this superficially describes all 2iV-photon scattering interactions in D di- 
mensions, we should note that in some cases the T-function may be undefined, i.e. for 
large even D and small N. This is not a serious problem as in the cases of most inter- 
est (e.g. 4 photon scattering in D = 4) the result will be either finite or manageable 
via regularisation. 

As a particular example we can consider the case of four photon scattering (N = 2) 
in four dimensional space-time. In this case the constraint over the {n^} implies that 
the only two terms allowed will have n\ = 2 and rii — 1. Thus we immediately recover 
the electromagnetic field invariants 

(Tr F 2 ) 2 , and Tr F 4 , 

obtained in previous discussions of four photon scattering. The coefficients can be 
evaluated easily from the general expression, and in this case we obtain 
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which is precisely the low energy limit of the result obtained by Schwinger [pi |i~0"f. 
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2.2 Spinor Case 

For this case our starting point is the spin- 1/2 1-loop Lagrangian of Brown & Duff 
(Eq. |2|) where we see the spinor case introduces the extra trace over spinor indices, 
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Before proceeding to use the techniques developed in the scalar case this trace must 
be evaluated and in order to do so it is neccessary to consider the cases of odd and 
even dimensional space-times separately. 



2.2.1 Even D 

Evaluation of the spinor trace is essentially a purely mathematical exercise, and in 
order to simplify the process as much as possible we can perform the calculation in 
Euclidean rather than Minkowskian space. In this case the metric, g^, assumes the 
form of the identity matrix and for even dimensions the antisymmetric electromag- 
netic field tensor can be written in the following canonical form, 
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where Aj are the eigenvalues. This canonical form can be achieved by a similarity 
transformation, which will not effect the commutation relations between the Lorentz 
group generators, but corresponds to choosing some canonical basis for the 7-matrices. 
If we consider a coupling of the electromagnetic field to a spin-1/2 field via 
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then this will result in the invariant matrix a ■ F which we can represent as 
\oF = h lfll „F^ = t( l0ll F 01 + --- + 7D _ 2lD ^ 1 F D - 2D - 1 



Noting that (7^) 2 = 1, i.e the identity matrix, the defining anticommutation relations 
imply 
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Consequently we can perform an expansion of the following exponential 
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which, on consideration of the canonical representation for the tensor F^ v ', gives 
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This result allows the exponential function of the full invariant matrix a ■ F, 
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where incidentally all the matrices in the exponent commute with each other, to be 
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Hence if we take the trace of this expression we obtain 
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where the factor of 1/2 which enters in the last line is due to the form of the canonical 
representation for F^ v . Substituting this evaluated form of the spinor trace into the 
effective Lagrangian results in 
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This is precisely the form required to use the techniques developed in the scalar 
case. We will proceed with this process later on, although now we should consider 
evaluation of the spinor trace for odd dimensional space-times. 
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2.2.2 Odd D 

In Euclidean space the canonical representation for the electromagnetic field tensor 
is given by the even dimensional case augmented by a zero eigenvalue Ao = 0. This 
is a consequence of the form of the characteristic identity for an odd- dimensional 
antisymmetric matrix which implies the existence of at least one vanishing eigenvalue. 
If we label the eigenvalues A = 0, Ai, A 2 , . . . , Am/2] then the analysis can proceed 
in an analogous manner to the even dimensional case. i.e. we obtain 
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where the extra eigenvalue has been included for clarity. Our labelling of the eigenval- 
ues now means that this expression is in exactly the same form as the even dimensional 
case and thus the interaction Lagrangian 
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applies in a space-time of arbitrary dimension. 

We may wonder at this stage as to the nature of any parity violating terms which 
may appear in odd dimensional space-times. These terms actually arise via the free 
field Lagrangian and must be considered separately, as we shall subsequently do. 



2.2.3 Evaluation of the Effective Interaction Lagrangian 

Having performed our manipulation of the Brown & Duff spin-1/2 1-loop Lagrangian 
into the form given by Eq. [7] we can follow our discussion of the scalar case, extracting 
a given order of the interaction (i.e. 2iV"-photon scattering), and expressing it in the 
following manifestly covariant form 
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The technique of Delbourgo & Matsuki jTT] can now be used to determine the coef- 
ficients a m in an analogous manner to the calculation in the scalar case. 
Following this procedure we have 
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The integration in this case makes use of the following identities 
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where, again, the symbol B n represents the Bernoulli numbers. We obtain in this 

case 
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Hence, for a 2iV-photon process mediated by a virtual spin-1/2 fermion loop, the 
interaction Lagrangian is given by 
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The final integration can be performed as in the scalar case giving 
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The expression comes of course with the proviso that the T-f unction is finite or can 
at least be handled via regularisation. 

We can obtain a compact form for the effective Lagrangian by noting the following 
correspondance between the Bernoulli numbers and the coefficients of the Bernoulli 
polynomial jTJJ], 
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With this substitution the coefficients are given by 

(_l)™+i2 4m (1 
a ™ ~ 4m(2m)! 2m \2 
and we can write the effective interaction Lagrangian in the following form 
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Hence the full effective Lagrangian, obtained by adding in the free field Lagrangian, 
is given in almost all dimensions by 
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The qualification is placed on this statement due to the existence of certain parity 
violating terms present in odd dimensional space-times which are expected due to the 
presence of fermion mass terms in the classical action. Delbourgo & Waites [15] have 
shown that these terms have the form 

Due to the order of this expression in the electromagnetic field tensor F, it is apparent 
that this term can only relate to a 2iV-photon process such that 

D = 2(2N) - 1. 

The relationship implies that in only one dimension D, for a given 2iV-photon inter- 
action, will the parity violating term emerge, i.e. for a 2 photon process D = 3, for 
a 4 photon process D = 7, etc. Interestingly, these are the dimensions for which a 
charge conjugation operator exists. Delbourgo & Waites have determined the general 
form of the coefficient C which, for an interaction of 2N soft photons, is given by 
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Hence, for such a process in a space-time of dimension D = 2(2N) — 1 (N integral of 
course), the effective Lagrangian is given by 
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In practice the only term of major importance is likely to be that of two photon vac- 
uum polarisation in D = 3, the extra invariant in this case being the Chern-Simons 
term. It is however useful, for reasons of generality, to have this all encompassing 
result. The form of such topological mass terms has been considered by Coleman & 
Hill H in QED 3 and by Reuter & Dittrich [TJ, and Delbourgo & Waites |TJ in 



arbitrary dimensional space-time. It was shown by Coleman & Hill in QED3, and 
in arbitrary dimensions by Delbourgo Sz Waites, that the particular parity violating 
term considered is indeed the only possible form allowed for reasons of gauge invari- 
ance. This results in the strict correspondence between the number of photons in the 
interaction and the dimensionality of space-time. 

Returning to a more general discussion we can observe that the interaction La- 
grangian is very similar to the expression obtained in the scalar case, with the only 
significant discrepancies being the usual spinor factor 2^/ 2 ], the different field invari- 
ant coefficients, and the possibility of parity violating terms. The introduction of the 
Bernoulli polynomial coefficients to describe the coefficients of the field invariants has 
exhibited the following correspondance: 

-£>2m (0) for spin-0 propagator particles, 

2 2m i?2m (-^) for spin-1/2 propagator particles. 



2 

It is perhaps a rather large leap of faith to suggest that the Bernoulli polynomial 
parameter is related to the spin of the virtual propagator particles, but at least this 
does introduce a very compact and illuminating representation. 

As a particular example we can consider the case of four photon scattering (N = 2) 
in arbitrary dimensional space-time. In this case the constraint over the {n^} implies 
that the only two terms allowed will have n\ — 2 and n% — 1. Thus we immediately 
recover the standard electromagnetic field invariants 

(Tr F 2 ) 2 , and Tr F 4 . 

The coefficients can be evaluated easily from the general expression, and in the case 
of four photon scattering we obtain 
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Using a technique due to Ackhiezer & Berestetskii |L8[ this result was verified from 
first principles by evaluating part of the Feynman amplitude for four photon scattering 
in arbitrary dimensional space-time. If we take the four dimensional case we obtain 
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which is the standard result obtained by Karplus & Neuman || and Schwinger [[|. 

This concludes the discussion of 2iV-photon scattering and to summarise we have 
obtained expressions (Eq. for the effective Lagrangian of the electromagnetic 

field which can describe photon-photon scattering to any order in arbitrary dimen- 
sions, where the virtual propagator particles may be either spin-0 or spin-1/2 particles. 
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3 2N+l-Photon Case 

The case of interactions between an odd number of photons, although not perhaps 
describing pure scattering, is worthy of comment to complete the discussion. If we 
consider the case of an odd number of photons interacting via a virtual fermion loop 
then the case of even D is trivial and results in a vanishing interaction Lagrangian as 
a consequence of Furry's Theorem. 

For odd dimensions the situation is not quite so simple as Furry's theorem no 
longer applies when D = 4k + 1 and there exists the possibility of having parity 
violating invariants similar in form to those which emerged in 2iV-photon interactions. 



Delbourgo & Waites |15[] have shown that these terms are given by 



2{2N + 1)!(4tt) 2 ^ 



for a 2N + 1-photon process (N integral). For a given 2N + 1-photon interaction this 
term is only present in a space-time with dimension D given by 

D = 2{2N) + 1, 

which is the set of dimensionalities for which the charge conjugation operator does 
not exist. Consequently, we can write the effective electromagnetic field Lagrangian 
in this case as 
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with the term of primary interest corresponding to a 3 photon interaction in 5D when 
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As mentioned earlier Coleman & Hill [16], and Delbourgo & Waites [0, have shown 
that these are the only possible terms which may arise, due to gauge invariance 
considerations, and thus interactions of an odd number of photons in odd dimensional 
space-times only have a non-zero effective interaction Lagrangian in the specific cases 
described above. 



4 Relations Between Invariants 

When considering effective Lagrangians in arbitrary space-time dimensions it turns 
out that the electromagnetic field invariants can suffer some linear dependence in 
certain dimensionalities. This property is most easily illustrated in terms of a partic- 
ular representation for the characteristic identities of the electromagnetic field tensor. 
This representation involves writing the characteristic equation in terms of polyno- 
mials over traces of matrix powers. A Maple routine has previously been developed 
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to obtain these results for any given dimension ||19|| . Some examples are listed below. 




D = A F 4 - \ T 2 F 2 + (§T 2 - |T 4 ) = 

D = 5 F 5 - ± T 2 F 3 + (|T| - ±T 4 ) F = 

D = 6 F 6 — | T 2 F 4 + (|Tl - |T 4 ) F 2 - ^T 3 + §T 2 T 4 - ±T 6 = 0, 



where Ti = Tr(F*). 

If we take as an example a four photon scattering interaction, the two electro- 
magentic field invariants, Tr F 4 and (Tr F 2 ) 2 , which are independent in D = 4, are 
linearly dependent in D = 2, 3. As can be seen from the characteristic identities 
above, in these cases 



Consequently the four-photon effective interaction Lagrangian simplifies in the case 
of D = 2 to 



Similarly higher dimensional cases can have any linear dependence in the invari- 
ants removed by using the relevant characteristic identity, the general condition for 
linear dependence being simply that the number of photons in the interactions is 
greater than the space-time dimension. This is precisely the situation in which we 
would expect some linear depenence of the photon momenta due to the limited num- 
ber of degrees of freedom. Thus it appears that the correspondence between F^ v and 
the photon momenta k^, via F^ v = {(k^A^lk) — k v ' A^ik)), causes any linear depen- 
dence in the photon momenta to manifest itself as a linear dependence in the field 
invariants. 

It is clear that at least in the simpler cases this linear dependence can be used to 
simplify the form of the effective Lagrangian. It would be interesting to know if the 
corresponding linear dependence in the photon momenta could be used to simplify 
calculation of the corresponding Feynman amplitude. 



(Tr F 2 ) 2 



2Tr F 4 . 




and in the case of D = 3 to 




5 Conclusions 



The discussion in this paper has focussed on low-energy photon interactions mediated 
by scalar and spinor propagators. We may also consider the possibility of interactions 
mediated by other forms of propagator particles. Photon-photon scattering mediated 
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by electroweak propagators such as W + and W particles has been considered by 



Okhlopkova [20|, and the gauge invariance of the amplitude in the electroweak the- 



ory has been established by Boudjema [21]. However, there has apparently been no 



discussion of the invariant form of the Lagrangian in the literature. We noted earlier 
a possible correspondence in the Lagrangians between the Bernoulli polynomial pa- 
rameter and the spin of the propagator particles. This relationship could be further 
tested by considering the form of the effective interaction Lagrangian for a photon 
scattering process mediated by a vector boson loop, i.e. W ± spin-1 particles. Taking 
the correspondence literally we might expect the effective Lagrangian to have a sim- 
ilar form to the expressions obtained for the scalar and spinor cases, except that the 
invariants would have coefficients involving, in some manner, the constants 

B2m (1) • 

At this stage, this is obviously only pure supposition but continuing in the same 
vein we may be able to consider a slightly more believable generalisation. The effective 
Lagrangians obtained by Brown & Duff || for the scalar and spinor cases indicated 
that essentially the only difference in the spinor case was the addition of an extra trace 
over spinor indices. Consequently we might predict that in the vector case similar 
expressions may result with the trace over spinor indices replaced by a trace of the 
form 



Tr exp [esS ■ F] , 

where £ M are the matrices equivalent to the Dirac matrices but for spin-1 particles. 
Realistically there may be Fadeev-Popov ghost contributions to take into account 
here, but the form of the 1-loop Lagrangians is certainly suggestive of the effective 
Lagrangian for photon scattering mediated by any electroweak propagator particle 
having a similar form to the expressions obtained in this paper. 
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